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Abstract 

We consider the problem of hypotheses testing with the basic sim- 
ple hypothesis: observed sequence of points corresponds to stationary 
Poisson process with known intensity against a composite one-sided 
parametric alternative that this is a self-correcting point process. The 
underlying family of measures is locally asymptotically quadratic and 
we describe the behavior of score function, likelihood ratio and Wald 
tests in the asymptotics of large samples. The results of numerical 
simulations are presented. 

AMS 1991 Classification: 62M05. 

Key words: Poisson process, self-correcting process, hypotheses testing, Wald 
test, likelihood ratio test, score function test. 



1 Introduction 

The model of self- correcting point process was proposed in 1972 by Isham 
and Wescott jTO] to describe a stationary sequence of events {ti, ^2, • • •} which 
does not have the property of Poisson process of independence of increments 
on the disjoint intervals. To introduce this processes we denote by X = 
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{^t, t >0} the counting process, i.e., Xt is equal to the number of events 
on the time interval [0,t]. Recall that for a stationary Poisson process with 
a constant intensity S > the increments of X on disjoint intervals are 
independent and distributed according to Poisson law 

P {Xt - Xs = k} = e-^^'-'\ 0<s<t, A; = 0,1,.... 

kl 

Particularly, 

P{Xt+^t~Xt>0} = Sdt (1 + 0(1)). 
For self-correcting point process we have 

P{Xt+^t-Xt>0\J^t} = S{t,Xt) dt (l + o(l)) 

where J-'t is the cr-field generated by {Xs, < s < t} and the intensity func- 
tion 

S{t,Xt) = aip{at-Xt), t>0. 
Here a > and the function ip {■) satisfies the following conditions: 

1. < ip {x) < oo for any a; G R, 

2. there exists a positive constant c such that (x) > c for any x > 0, 

3. lim ip {x) > 1, and lim ip (x) < 1. 

Self-correcting processes are called as well stress-release processes (see ji], 
p. 239). This class of processes is widely used as a good mathematical model 
for non-poissonian sequences of events. This model was found especially 
attractive in the description of earthquakes (see Ogata and Vere- Jones 
Lu at al mi). 

Example 1. Let 

S{t,Xt)=e^v{a + l3{t- QXt)} 

where /? > 0, ^ > 0. It is easy to see that the conditions 1-3 are fulfilled and 
the point process with such intensity function is self-correcting. 

This model was studied by many authors (see the references in 0]). Par- 
ticularly it was shown that under mild conditions there exists an invariant 
measure ix and the law of large numbers (LLN) 

{St -Xt)dt^ j h {y) IX (dy) (1) 

2 



is valid (see Vere- Jones and Ogata [17j, Hayashi |^9J), Zheng j21j. Here h{-) 
is a continuous, integrable (w.r.t. /i) function and 5 > is the rate of the 
point process. For the model of Example 1 we have the LLN if p > and 



As the self- correcting model is an alternative for the stationary Poisson 
process, it is natural and important to test these two hypotheses by the 
observations {^1,^2? • • •} on the time interval [0,T], i.e., to test 



Remind that the likelihood ratio in this problem has the following form 



where Xt_ is the limit from the left of Xt at the point t y^. Therefore, if 
the function atp {■) /S is separated from 1 then the second integral in this 
representation tends to infinity and there are many consistent tests. Hence 
it is more interesting to compare tests in the situations when the alternatives 
are contigous, i.e. the corresponding sequence of measures are contigous. 
This corresponds well to Pitman's approach in hypotheses testing [19j. We 
can have such situations if ip {■) = S + o(l) with special rates o(l). In 
this work we consider one of such models defined by the intensity function 
5* (t, Xt) = Sip {'d {St — Xt)) where is a small parameter and ip (0) = 1. We 
suppose that the function ip (■) is smooth and we can write 



It is easy to see that the rate = — ^ under hypothesis S {t, Xt) = is 
~ because 



(3>0. 



S (t, Xt) = S versus S (t, Xt) = aip {at — Xt) . 
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where Wt {s) = {ST)'^^^ {STs - Xts) ^W{s), and {W (s) , < s < 1} is 
Wiener process. Note that we put a = S, otherwise 



{at-Xtf dt 



^ (0)^ 4 



^(0) 
6 



T [{a-S)vT+^STWT{v)y dv 



^l{a-Sf (1 + 0(1)) 



Therefore, if a 7^ 5, then we have to take = uT~^^'^ and to test the 
simple hypothesis : m = against Mi : m > 0. In this case the family 
of measures is LAN and the usual construction provides us asymptotically 
uniformly most powerful test (see, e.g., Roussas [2III)- Note that according 
to for any fixed alternative "i? > we have the convergence 



\ j\st - Xtf dt I y'fiidy) 



T 

which, of course, requires another normalization. 

Therefore we consider the problem of hypotheses testing when under hy- 
pothesis Mo the intensity function is a known constant S > (Poisson pro- 
cess) and the alternative ^ is one-sided composite: self-correcting process 
with intensity function S (t, Xt) = Sip i'&T {St — Xt)), where for convenience 
of notation we put {}t = u/S'ip{Q)T (we suppose that ip {Qi) > 0). In this 
case the corresponding likelihood ratio Zt (u) converges to the limit process 

Z{u)=exp\-u [ W{s) dW{s)~\ f W{sfds 



i.e., the family of measures is locally asymptotically quadratic [12] . We study 
three tests: score function test, likelihood ratio test, Wald test and compare 
their power functions with the power function of the Neyman- Pearson test. 
Note that we calculate all limits under hypothesis (Poisson process) and we 
obtain the limit distributions of the underlying statistics under alternative 
(self-correcting process) with the help of Le Cam's Third Lemma. Therefore 
we do not use directly the conditions 1-3 given above. 
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The similar limit likelihood ratio process arises in the problem of hypothe- 
ses testing u = against m > for the time series 

Xj = (l-^^ Xj^i + ej, j = l,...,n^ oo, 

where ej are i.i.d. random variables, 'Eej = 0, Ee^ = a^. The asymptotic 
properties of tests are described under hypothesis and alternatives by Chan 
and Wei |2I and Phillips [18 . Particularly, the limits of the power functions 
are given with the help of Ornstein-Uhlenbeck process 

dYs = -u Ys ds + dWs, = 0, < s < 1. 

Then Swensen [22] compared these limit powers. 

For the model of Example 1 the power functions (for local alternatives) 
was studied by Ogata and Vere- Jones |l22j and by Luschgy [I^, (TH]. The 
limit likelihood ratio and tests are similar to that of the mentioned above 
time series problem. Remind as well that Feigin [6] noted that the same 
limit likelihood ratio arises in the problem of testing the simple hypothesis 
M = against one-sided alternative m > by observations 

dXt = -- Xt dt + dWt, Xo = 0, < t < T ^ cx). 

In our case we obtain similar limit expressions for the likelihood ratio and 
power functions and compare the errors of tests. The analytical considera- 
tions give us an asymptotic (for large values of u) ordering of the tests. The 
numerical simulations of the tests show that for the small values of e and for 
the moderate values of u the power functions of the likelihood ratio and Wald 
tests are indistinguishable (from the point of view of numerical simulations) 
of the Neyman-Pearson envelope. This interesting property was noticed (for 
e = 0.05) by Eliott at al. j5] on the base of 2 ■ 10^ simulations. In our work 
we obtain similar result having 10^ simulations and we observe for the larger 
values of e that the asymptotic ordering of the tests holds already for the 
moderate values of u. 

A similar problem of hypotheses testing in the situation, when the alter- 
native process is self-exciting ^ was considered in jS]. 

2 Score Function Test 

We observe a trajectory X'^ = {Xt, < t < T} of a point process of intensity 
function S {■,Xt) and consider the problem of testing the simple hypothesis 
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against close one sided composite alternative 

: S (t, Xt) = S^, (2) 

^i: S{t,Xt)^S,ilj{i9T[Sj-Xt]), i9t>0, (3) 

where is a small parameter, the value and the function ip {•) are known. 
The problem is regular in the following sense. 

Condition A. The function -0 (x) , a; G IR is positive, continuously diSer- 
entiable at the point x — 0, t/j (0) — 1 and ip (0) > 0. 

The rate of convergence ~^ is chosen such that the likelihood ratio 
L is asymptotically non degenerate. In the case ip {0) < we need 

to change just one sign in the test. This leads us to the reparametrization 

?/ 

= -■ , u>0 

and to the corresponding hypotheses testing problem 

J%: u^O, (4) 

JTi : u>0. (5) 

Therefore, we observe a Poisson process of intensity -S* under hypothesis 
and the point process under alternative has intensity function 

S {t, Xt) = 5, + I {SJ -Xt)+o . 

Let us fix £ e (0, 1) and denote by Jt^ the class of test functions (f)T (X^) 
of asymptotic size e, i.e., for 0r e ^ we have 

lim Eo dr (X^) = s. 

T-^oo ^ ' 

As usual, (-^^) is the probability to accept the hypothesis M\ having 
observations X^ . The corresponding power function is 

{u, 0^) = E„ 4>T {X^) , > 0. 

Let us introduce the statistic 

At (X^) =^ r {SJ - Xt.) [dXt - dt] 
^* ^ Jo 



Xx — (yXx — SJT)"^ 
2S?f ■ 



(6) 
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This equality follows from the elementary representation (see, e.g. 
Lemma 4.2.1) for the centered Poisson process tt^ — — t 



TT. 



T 



Tit- dvTt + TTr + S^T. 



which obviously is equivalent to 



- I nt- dnt - 



2T 

Define as well two random variables 

W (s) dW (s) , 



W(sYds, 



AiW) = l{l-W{lf) = -J^' 

where {W (s) ,0 < s < 1} is standard Wiener process. 

Remind that the likelihood ratio in this problem has the following form 

m (7 = 5*^(0)) 

L (:^.^'^) = exp (S, t - Xt-)^ [dXt - dt] 

i;(^-^iS.t-Xt)^-l-ln^(^-^iS.t-X,)^ 5,dt|. (7) 

Therefore the direct differentiation w.r.t. u at the point u = gives us 
the introduced above statistic 



_d_ 

d u 



InL 



u 

if-' 



X' 



At (X^) . 



u=0 



Below we denote 



1 - zt. 



and 



h{u) 



2u 



l-e 



-2u ' 



where is 1 — a quantile of standard Gaussian law, i.e., P (C > Za) = a, for 
C~Ar(0,l). 

We have the following result. 

Theorem 1 Let the Condition A be fulfilled, then the score function test 
belongs to the class and for any u^, > its power function 



(9) 



7 



Proof. Under hypothesis the value is a poissonian random variable 
with parameter S'^.T. Therefore we have immediately 

and At (X^) ^ A (ly) as T ^ oo. Hence 

Po {At (X^) > aj P |a W > = P {iCI < = ^. 

This provides </)t G 

To study the power Pt{u*,(Pt) "we would like to use the Third Le Cam 
Lemma ^21; [21] ■ Therefore we need first to show the joint weak convergence 



Co (At, It (u)) =^ C ( A{W), uA{W) - y J(W^)) (10) 



where h (u) = InL (:^,^^)- 
To verify (fTUjl we denote 



where 



and show that 



(u) =uAr (X^) - ^ Jt (X^) 



Jt (x^) =^,[ (s.t - X,f dt 



Co {It («)) =^ ^ (^A(W^) - y J(W^)) • (11) 

Then fllOj) will follow from the convergence 

/t («t) - (mt) (12) 
for any bounded sequence ut- 

Lemma 1 

Co {At (X^) , Jt (X^) ] ^ j\ {s) dW {s) , W (.)' ds^ . (13) 
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Proof. Let us put Wt (s) = (5,T)"^/V,t, s e [0, 1]. Then 

EoWr (s) = 0, Eo [Wt {si) Wt (ss)] = min (si, Ss) 

and we have 

Jt {X^) = ^ £ vr^ dt = j'^ Wt {sf ds. 

Using the standard arguments we verify (well-known fact) that for any col- 
lection {si, . . . , Sk] we have the weak convergence (as T oo) of the vectors 

{wT{s{),...,WT{srS) =^ {w{s^),...,W{sk)). 
Moreover the following estimate holds 

{Eo\WT{s,f-WT{s2f\Y< 

< Eo \Wt (si) - Wt (s2)|' Eq \Wt (si) + Wt (s^)]^ < 4 |s2 - Si| . 

Hence (see Gikhman and Skorohod [Zj, Section IX. 7) we have the conver- 
gence (in distribution) of integrals 

Wt(s? ds^ [ W(s? ds 



and 

A. (XT ^ i^iMl£ ,1 , „ ^ . _ I (,) (,) . 

It is easy to see that we have the same time the joint convergence too because 
from the given above proof it follows that for any Ai, A2 

XiWt {i f + A2 / Wt {sf ds =^ XiW {if + X2 [ W {sf ds. 
Jo Jo 

Therefore the Lemma ^ is proved. 

Our goal now is to establish a slightly more strong than (fT^ relation 

Ij, (ut) = UT At (X^) (1 + o (1)) ['wt {sf ds {1 + o (1)) (14) 

Jo 

where o (1) ^ for any sequence Ut G Uy with Dt = {u : < u < ^•^*^ }. 



We can write 



u Wt (I) 



— In'^ 



-uWt (^) 



T 



« 5l,T - — S2,T 



dTTt 

-uWt (^) 



dt 



with obvious notation. Remind that u > 0. Using Lenglart inequahty we 
obtain for the first term 



Po{\Si,t\ >a}<- 
a 



Po 



' r ,/s:jT , , / -uWt {s) 
Wt (s) + In ^/^ ' ^ ^ 



u 



ds > b 



for any a > and 6 > 0. Now expanding the functions {■) we obtain 



= 1 



uWt (s) ^ /^mH^t (s) 



where it < i*. Introduce the set 

Ct^Iou: sup |Wt(s)| < -0(0) Vhi^l 

I 0<s<l J 

and note that for a; e Ct we have the estimate 

u\Wt(s)\ ^ 1 
sup sup -. < 



ueVr 0<s<l 

Hence for all m G Ut on this set wc can write 

-uWt (s 



sup 

0<s<l 



V'(O)-^ 



< sup 

|t)|<(lnr)-l/2 



V'(o)-V'(^) 



= /ir ^ 



as T — > oo because the derivative is continuous at the point v — 0. 
Let us denote Ug — 1^77^^- Using the expansion of the logarithm 

In i-us)) = In (1 - u^i^ i-u,)) = --^4t^- 
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we obtain the following estimate 

»i 







Wt (s) + In ^ {-Us 



u 



+ Po 



ds > S> < Po {C^} 

i) {-Us) 



Wt {sY 1 1 

\ ^(0) [l-Usi^{-Us) 



ds > 6, Ct > . 



Remind that Wt {s) is martingale, hence by Doob inequality we have 

Po{C^} < Po {\Wt (1)1 > ^ (0) VwT] < J 
For the second probability after elementary estimates we obtain 

2 

ds > 6, Cr > < 



1p {-Us) 



Po < / Wt {s)^ 1 - — 

lo \^ ^(0) {l-u^i;{-Us) 



< Po <^ C / Wt {s)" ds + 



^ - ''-2 ' M >a <£ if^2 







InT 



2b 



1 



T 



InT 



with some constant C > 0. Recall that by Tchebyshev inequality 







WT{s)'ds>A\< ^. 



Therefore, if we take h = then for any a > 

Po{|5i,t| >a}^0 

as T ^ CX3. 

The similar arguments allow to prove the convergence 

Po{IVI > a} ^0 

too. 

Therefore, the likelihood ratio Zt {u) = L , u > is (under 

hypothesis Mq) locally asymptotically quadratic (LAQ) [T2j, because 

Zt {u) ^ Z {u) = exp !^-u W {s) dW {s) - ^ W {s)^ ds^ (15) 
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Moreover, we have the convergence {ut) — It (wt) ~^ for any bounded 
sequence of € Ut- Note that the random function Z {u) is the hkehhood 
ratio in the hypotheses testing problem 

JTo: u = 0, 

JTi: u>0, 

by observations of Ornstein-Uhlenbeck process 

dY (s) = -uY (s) ds + dW (s) , Y{0) = 0, 0<s<l (16) 

under hypothesis u = 0. 

This hmit for the hkehhood ratio under alternative can be obtained di- 
rectly as follows. Let us denote 

Then using the representation 

Xt = S [ ^{^T[S,r-Xr]) dr + Mt 
Jo 

where Mt is local martingale and expansion of the function ip (■) at the vicin- 
ity of we obtain the equation 

Yt (s) = ~u I t^Yr (f ) dv + Vt (s) , Yt {0) = 0, < s < 1 
Jo ip (0) 

where Vr (s) is local martingale and gy = ■^^^^^^=Yt iy) 0. The cen- 
tral limit theorem for local martingales provides the convergence Vt (s) =^ 
W (s). Hence the process (fTBj) is the hmit (in distribution) of Yt (s). More- 
over from ((Tj) we have 

AT(xn=^ + i^^^i^^. 
^ ' 2^/S7T 2 2 

This limit of the statistic At (-^^) follows from the Third Le Cam Lemma 
as well. Particularly, for any continuous bounded function H (■) 

E„/7 (At (X^)) = Eo [Zt (n) H (At (X^))] 

Eo [Z {u) H (A (W))] = EuH (A (F)) , 
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where 



A{Y) 



Y(s) dY(s) 



l-Y(l) 



2m* 



Hence under alternative {-^t = u^/'-yT) we have the convergence 
(5t{u,At)^^u,{\Y{1)\ <zi^\=v[\W{l)\ <z, 
because 



r(i) 

This proves Q. 



e-"(i-^) diy(s) -ATIO 



1 



-2m, 



Theorem n is asymptotic in nature, and it is interesting to see the powers 
of the score function test for the moderate values of T and especially to com- 
pare them with the hmit power functions. This can be done using numerical 
simulations. 

We consider the model of Example 1 with S*.,, = 1 and il){t) = e*. This 
yields the intensity function 



S{u,t,Xt)=exp(^^ [t-Xtij, 



u>0, 0<t<T. 

In Figure 1 we represent the power function of the score function test 0^ 
of asymptotic size 0.05 given by 

{u, = Pm {At (X^) > ao.05} , < M < 20, 

for T = 100, 300 and 1000, as well as the limiting power function f3*{-) given 
by the formula 0. 

1 

0.9 
0.8 
0.7 
0.6 



t o.s 



a. 

0.4 



T= 100 
T= 300 
T=1000 
Limiting 




Fig. 1: Power of the score function test 
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The function 0^) is estimated in the following way. We simulate (for 
each value of m) M = 10^ trajectories Xj , j = 1, . . . , M of self-correcting pro- 
cess of intensity 5* {u,t,Xt) and calculate Aj = At{XJ). Then we calculate 
the empirical frequency of accepting the alternative hypothesis 

1 

— y^i,^ , ^ dju,(t)*). 

Note that for T = 1000 the limiting power function is practically attained. 
Note also that for T = 100 the size of the test is 0.079 which explains the 
position of the corresponding curve. 

Remind that score-function test is locally optimal p. 

3 The Likelihood Ratio Test and the Wald 
Test 

Let us study two other well-known tests: the likelihood ratio test (px based 
on the maximum of the likelihood ratio function and the Wald test 0t based 
on the MLE dj,. 

Remind that the log-likelihood ratio formula is 

In L {d, X^) = [ In ^ (t9 (Sj ~ Xt-)) [dXt - dt] 
Jo 

- [ [i^i'd {SJ - Xt^)) -l-ln^{^ {SJ - Xt_))] S, dt 
Jo 

and the likelihood ratio test is based on the statistic 

6t {X^) = supL , 

where 6 is the set of values of under alternative. The test is given by the 
decision function 

where the threshold bs is chosen from the condition (pT G J(fs- 

Note that 6t {X^) = L {§t^X^^ as well, where "dx is the maximum 
likelihood estimator of the parameter d. 

The reparametrization -d = "Ot = uj^T reduces the problem ©-© to 
(^-(jni) and we have to precise the region of local alternatives. In the tradi- 
tional approach of locally asymptotically uniformly most powerful tests [201 
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(regular case) to check the optimahty of a test 0^ we compare the power 
function Pt {u, (pj) with the power function of the Neyman-Pearson test on 
the compacts < m < for any > 0. For these values of u the alter- 
natives are always contigous. To consider the similar class of alternatives in 
our case is not reasonable because the constant 6^ became dependent of K. 
Indeed if we take the test function 

0^(X^)=lr Zt{u)=l(^,X^ 

i sup ZT{u)>b, } \11 
Lo<ti<K J 

then the condition G imphes 6^ = 6^ (K). Therefore we suppose that 

Finally, we have the following hypotheses testing problem 

J^o: u = 0, (17) 

^i: u = u^eVT (18) 

Therefore, to study 

6r {X'^) = 1 

sup ZT{u)>b^ 



we need to describe the asymptotics of its errors under hypothesis Jifo and 
;rnativ( 
Below 



alternatives Jifi with {} = E Vt- 



A(W)- ^<^''> 



^/2J(W)' 

Theorem 2 Let us suppose that condition A is fulfilled and the value is 
solution of the equation 

P {k{W) > be) = e. (19) 

Then the test (j)T with = e^^ belongs to J^s OL^d its power function converges 
to the following limit 

f3 (m„ 0t) /3 (u,) = P {A(KJ > be} , 

where 



•/iWZ) V8 j(y'„.) 

and Yu, = {Yu, (s) , < s < 1} is Ornstein- Uhlenbeck process pfjj) with u 
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Proof. The log-likelihood process It {u) = In {u) admits (under hypoth- 
esis J^) the representation ([11)) 



It (u) =uAt (X^) (1 + 5i,t) " y Jt {X^) (1 + hr) (20) 



where Si^T uniformly on n G Vt- Hence 
At {X^Y = sup hiu) = 

u£Vt 



2J{W) 



and we have 

Eo0T {X^) = Po ( sup It (u) > bl] P (A(H^) > 6,) = e. 
Let us fix an alternative u = u^,. We have the convergence 



Co {At (X^) , It («*)} ^ /: |a(W), A(W) - ^ J(W)| . (21) 

The convergence ()21|) allows us to apply Third Le Cam's Lemma as fol- 
lows: for any bounded continuous function H {■) 

E^,H {At (X^)) = Eo [Zt (m,) H {At (X^))] 

Eo [Z (u,) H (A (W))] = Eu,H (A (F,J) . 

Hence 

f3 {u,, 0t) = Pn. ( sup It (u) > bA {A (F„J > 6J . 

This completes the proof of the theorem |21 

Let us note, that the threshold bs is given implicitly as the solution of 
the equation (fT^. In the following table we give some values of b^ obtained 
using numerical simulations. 



e 


0.01 


0.02 


0.03 


0.04 


0.05 


0.1 


be 


1.814 


1.636 


1.524 


1.440 


1.373 


1.144 



These thresholds are obtained by simulating M = 10'' trajectories on 
[0, 1] of a standard Wiener process, calculating for each of them the quantity 
A(iy) and taking (1 — £:)M-th greatest between them. 
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The next test usually studied in such hypotheses testing problems is the 
Wald test 

where -^t is the maximum likelihood estimator of 
Below 

Theorem 3 Let us suppose that condition A is fulfilled and the value is 
solution of the equation 

P iT{W) > c) = e. (22) 

Then the test (pT belongs to ^ and its power function for any alternative 
converges to the following limit 

(3 (w*, 0t) /3 K) = p {r(Kj > cj , 

where 



3(Yu.) J(K,) 

and Yu^ is the same as in Theorem\^ 

Proof. The proof follows immediately from the representation ()2U|). because 

Pf ^ i-iT^T >Ce\= Pf ^ I sup [u] < sup Zj, (u) \ > 



— ^ Po <^ sup Z {u) < sup Z{u)\ {T{W) > c J = e 
and (under alternative u = u^) 

PiP {iT^T >Ce]= PiP ( sup Z^ {U) < sup Z^ {U) \ 

p„, I sup z {u) < sup z (n) I = p {r(y;j > cj = (u,) 



As above, the threshold q is given implicitly as the solution of the equa- 
tion (j22I)- In the following table we give some values of q obtained using 
numerical simulations. 
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e 


0.01 


0.02 


0.03 


0.04 


0.05 


0.1 




13.692 


11.224 


9.803 


8.806 


8.042 


5.719 



These thresholds are obtained by simulating M = 10 trajectories on 
[0, 1] of a standard Wiener process, calculating for each of them the quantity 
T{W) and taking (1 — e)M-th greatest between them. 



4 Comparison of the Tests 

Remind that all these three tests 0^, (pT and 0t in regular (LAN) case are 
asymptotically equivalent to the Neyman- Pearson test 0° (with known al- 
ternative u) and hence are asymptotically uniformly most powerful. In our 
singular situation all of them have different asymptotic behavior and there- 
fore it is interesting to compare their limit power functions 

{u) = P„ {A (K) > a J , (3 (u) = P„ <; ^ffiL > b. 



2J(y;) 



,2 



of course, under condition that all of them belong to /C^. Our goal is to 
compare these quantities for the large values of u. 

We have to study the distribution of the vector (A (Yu) , J (V^^)), where 

A(Y^) = - [ Yu{s)dY^{s), Hyu)= [ Y^{s)' ds, 
Jo Jo 

where Yu is solution of the equation 

dYuis) = -uYuis) ds + dW (s) , y„ (0) = 0, < s < 1. 

Let us introduce the stochastic process = -Ju Y^ (^) , < f < m (this 
transformation was introduced by Luschgy 16ij). Then we can write 

^Vv = -Vv dv + dwy, 1/0 = 0, < f < M, 

where = y/u W (^) is a Wiener process and 

A (K) = -U-' r dy, ^ ^, J (K) = u-^ f yl dv ^ \ 
Jo u Jo 
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in obvious notation. Further, the process i/y is ergodic with the density of 
the invariant law / (y) = / ^/tF. Hence J„ ^ oo and 

Note that the distribution of the process does not depend on u. 
The constant = {u) because it is defined by the equation 

uA{W)-—^W)>dA=e. 

For the large values of u this constant can be approximated as follows. We 
have (under hypothesis Ho) as w — > oo 

P„|i,A(H')-^J(H')>4(«)| = 



u 

1 



where the constant is defined by the last equality. For example, if we take 
e — 0, 05 then the numerical simulation gives us the value eo,o5 = 0, 056. 
Therefore d^ (u) — —0, 5 eg (1 + o (1)). If we suppose that e is small and 
try to solve the equation 

/ /j (x) dx^£ 
Jo 

where /j (x) is the density function of the integral J (W), then we can easily 
see that /j (0) = and all its derivatives /j^'' (0) = 0, A; = 1, 2, . . .. Hence 
to see an approximative solution we need to calculate the large deviation 
probability of the following form (below r = s/^/e^, E — Ce^^'^ ^ oo). 

Poje;^^ W^(s)'ds < l| =Po|^ W{rfdr<l^. 
Below we put d^ {u) — —0, 5 u^. 
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We have the relations 

P* (n) = P { A„ > n aj = P jy dw^ < J„ - n| , 

P{u)=p\^>-]=pI [ y,dw,<3u--3u 
13° {u) = P |a„ - ^ > 4} = P y\.dw, < + 
Therefore the large values of u (J„ ~ u/2) 

> 3u — be V 2J„ > 3u — O-s 

2 2 u 

and finally 

/?* (m) < ^ (m) < /3 (m) < /?° (m) . 

These inequalities are in accord with [22] • 

Note that for small values of e the constant is close to 0,5 (e.g. ao,o5 = 
0,498, ao,oi = 0,49992) and in this asymptotics the power of score-function 
test is 

f3*{u) = p\ [ y,dw, < (0,5-a,) M(l + o(l))j. 



Hence one can expect that in this case the score-function test has essentially 
smaller power than the others. 

Now let us turn to numerical simulations of the limiting power functions. 
We aim to obtain the limiting power functions of all the three tests, as well 
as the Neyman-Pearson envelope, for the moderate values of u {u < 15). 

Note that for the score function test P*{u) can be computed directly 
using 0. However the limiting power functions of the likelihood ratio and of 
the Wald tests are written as probabilities of some events related to Ornstein- 
Uhlenbeck process and can be obtained using numerical simulations. 

For the likelihood ratio test we have 

where 



Z (u) = exp |mA(W^) - y J(Vr)| . 
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So we simulate M = 10^ trajectories Wj = {Wj{s), < s < 1}, j = 
1, . . . , M of a standard Wiener process and calculate for each of them the 
quantities A^- = A{Wj), Jj = J(l^j), Aj = ^j/^j and (for each value of u) 

Zj (u) — exp jwAj — Then we calculate the empirical mean 

1 ^ 

For the Wald test we have similarly 
1 

where Tj = Aj/ a/SJj. 

Finally, in order to compute the Neyman-Pearson envelope, we first ap- 
proximate (for each value of u) the quantity = de{u) by the (1 — £)M-th 
greatest between the quantities \nZj{u), and then calculate 

1 ^ 

The results of these simulations for e — 0.05 are presented in Figure 2. 




u 



Fig. 2: Limiting powers for e = 0.05 

Let us note here that in this case the power functions of the likelihood 
ratio test and of the Wald test arc indistinguishable (from the point of view 
of numerical simulations) from the Neyman-Pearson envelope. This quite 
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surprising fact was already mentioned by Eliott at al. |5j, who showed the 
similar pictures having 2 ■ 10^ simulations. As we see from Figure 2, with 10^ 
simulations the curves are still indistinguishable. The situation is however 
different for bigger values of e. The results of simulations for e = 0.01, 0.05, 
0.25 and 0.5 are presented in Figure 3. 




One can note that for big values of e (e.g. e = 0.5) the powers became 
more distinguishable, and that the asymptotically established ordering of 
the tests holds already for these moderate values of u. Note also that for the 
small values of e (e.g. e = 0.01 and 0.05) the curve of score-function test is 
essentially lower as expected. 



5 Discussion 

Remark 1. Note that alternatives u = ut —>■ oo with i^uT ~^ ^'^^ local but 
not contigous. That means that the corresponding sequences of measures 
^P^"^^ , Pq^'' j , T ^ oo are not contigous. Particularly, the second integral in 
the likelihood ratio formula tends to infinity: 

f [ij {^u^ {SJ - Xi_)) - 1 - In {^^^ {SJ - Xi_))] dt ^ oo. 

In such situation the power function of any reasonable test tends to 1 and 
to compare tests we have to use, say, the large deviation principle. For 
example, the likelihood ratio test 0^ is consistent for the local far alternatives 
= /g^rp ,v G [u, V] where < u < V < oo. Indeed, under mild regularity 
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conditions we can write 



E,0; {X^) = Po ( sup L ( 



= Po { sup 

^u<v<V 

• 1 



= Po<; sup 

.i/<i;<y 

-1 



/ In V' (vWt (s) ) dWr (s) - 
Jo 

S,T I [V' {vWt (s)) - 1 - In {vWt {s))] ds 
Jo 

lmP{vWT{s))dWT{s)- 

- [ [ip {vWt (s)) - 1 - In {vWt (s))] ds 
Jo 



> In Ce > = 



— ^p| |nf^^ [^/^(wiy (s)) - 1 -lnV^(viy (s))]ds > o| = 1 

because the function g (y) = y — 1 — \ny > for y ^ 1 and g {y) — y — 1. 

Remark 2. Note, that we can construct asymptoticaUy uniformly most 
powerful test if we change the statement of the problem in the following way. 
Let us fix some D > and introduce the stopping time 

TD = inf |t : 2 {S,t - Xtf S, dt > L>'| . 

Then we consider the problem of testing hypotheses 

J^o : S {t, Xt) = S"*, 

Jifi : S{t,Xt) = S,ilj{i»D[S.t-Xt]), 



> 



ip{0)D 

by observations X'^^ = {Xt, < t < td} in the asymptotics D ^ oo. Now 
the likehhood ratio Z^^ (u) = L will be LAN: 



u 



C-A/'(o,i) 



and the test = ^{A,^(x-D)>2e} ^^^^^ 



A,^ = 1 - XtS) [dX, - SAt\ 
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is locally asymptotically uniformly most powerful. 

The proof follows from the central limit theorem for stochastic integrals 
and the standard arguments (for LAN famihes) . 

Remark 3. Note that these problems of hypotheses testing are similar to 
the corresponding problems of hypotheses testing for diffusion processes. In 
particular, let the observed process X'^ — {Xt, < i < T} be diffusion 

dXt = ijj i-^xXt) dt + adWt, Xq = 0, 0<t<T, 

where the function ^p (0)=0, is continuously diffcrcntiable at the point and 
ip (0) > 0. If we consider two hypotheses: -d — and -d > then the 
reparametrization 

ua 

Vt = -■ 

V'(O) T 

provides local contigous alternatives, i.e., the log- likelihood ratio in the prob- 
lem 

J^i : u>0. 

has the limit: 

In L f . , xA ^-u f'w (s) dW (s) - — f'w (sf ds. 

VV'(o)r J Jo ^ ^ 2 7o 

The score function test based on the statistic 

{X^)^-^j\,dXt, 

the likelihood ratio test and the Wald test have the same asymptotic prop- 
erties as those described in Theorems 1, 2 and 3 above. 

For example, if (x) = x, then we have the Wiener process (under hy- 
pothesis Mq) against ergodic Ornstein-Uhlenbeck process under alternative 

Remark 4. We supposed above that the derivative of the function [x] 
at the point x = is not equal to 0, but sometimes it can be interesting to 
study the score function and the likelihood ratio test in the situations when 
the first k — \ derivatives with k>2 are null. 

Let us consider a self-correcting process X^ = {Xt, 0<t<T} with inten- 
sity function S^ijj {Sj - Xt)) such that ^ (0) = 1 (0) = and (•) 7^ 
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{k = 2). In this case the modifications have to be the following. Suppose 
that t/j (0) > 0. To have LAQ family at the point -& — we chose the 
reparametrization i9 = i?„ 



which provides the limit 

lnL{^^,X^) / W{sfdW{s)-— / W{s)^ds 

Jo 2 Jo 

Then in the hypotheses testing problem 



Mi: u>0 



the score function test ip (^^) = (xT)>ce} based on the statistic 

{X^) = ^^^^ £ {SJ - X,f [dX, - SAt] . 
It is easy to see that under Mq 

A^(X^)^^-j[|V(.) d.. 
Hence to chose the threshold Cg we have to solve the following equation 



x'^—y>3c 



because (1) , 3 JqW (s) ds^ is Gaussian vector. 
The cases k > 2 can be treated in a similar way. 
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